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Abstract 

Chemotaxis along helical paths towards a target releasing a chemoattractant is found in sperm cells and 
many microorganisms. We discuss the stochastic differential geometry of the noisy helical swimming path 
of a chiral swimmer. A chiral swimmer equipped with a simple feedback system can navigate in a concen- 
tration gradient of chemoattractant. We derive an effective equation for the alignment of helical paths with 
a concentration gradient which is related to the alignment of a dipole in an external field. We discuss the 
chemotaxis index in the presence of fluctuations. 
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Biological microswimmers use flagellar propulsion or undulatory body movements to swim at 
low Reynolds numbers In addition to forward propulsion with translational velocity v, any 

chirality in swimming stroke results in a net angular velocity Q. Hence, such a swimmer moves 
along a helical path with curvature k = | 17 xv|/|v| 2 and torsion t = |fi-v|/JvP in the absence of 
fluctuations J3|l. Helical swimming paths have been observed for sperm cells HtslsLZQ, eukaryotic 
flagellates [8|,|90, marine zooplankton HE, and even large bacteria [| 1 2fl . A necessary condition 
for a pronounced helicity of the swimming path of a chiral swimmer is given by | f2 1 3> D rot where 
the rotational diffusion coefficient D rot ~ L~ 3 depends strongly on the size L of the swimmer 



[| 1 3L Il4fl . Thus there is a critical size for a chiral swimmer below which fluctuations diminish 
directional persistence and interfere with helical swimming. The bacterium E. coli for example 
is much smaller than the swimmers mentioned above and fluctuations dominate over an eventual 
chirality of swimming. Nevertheless, this bacterium can navigate in a concentration field of a 
chemoattractant by performing a biased random walk 111 311 . A larger swimmer moving along a 
helical path can exploit a fundamentally different chemotaxis strategy: It has been shown both 
experimentally [|il[aB&|9,12,15,16] and theoretically [17, 18] that such a chiral swimmer can 
navigate in a concentration gradient of chemoattractant by a simple feedback mechanism. Here we 
study the impact of fluctuations and show that sampling a concentration field along noisy helical 
paths is a robust strategy for chemotaxis in three dimensional space even in the presence of noise. 
The alignment of noisy helical paths with a concentration gradient is formally equivalent to the 
alignment of a polar molecule subject to rotational Brownian motion in an external electrical field. 

Stochastic differential geometry of noisy helical paths. The geometry of a swimming path r(i) 
is characterized by the tangent t = r/v, normal n = t/|t| and binormal b = t x n, where v — |r| 
igpeed and dots denote time derivatives. The time evolution of these vectors can be expressed as 



ft, t = vkii, n 



-vk t + vr b, b 



-vr n 



(1) 



where n(t) and r(t) are curvature and torsion of the swimming path r(t), respectively. For a noisy 
path, K,(t) and r(t) fluctuate around their mean values 



«(*) =«0 +£*(*), T{t)=T + £ T {t) 



(2) 



where £ K and £ T are stochastic processes with mean zero and respective power spectra S K , S T , as 
well as a cross power spectrum S K>T Q26Q . For simplicity, v(t) = v is assumed constant. The 
stochastic differential equations (|1I2|) involve multiplicative noise and should be interpreted in the 
Stratonovich sense if £ K or £ r is 5-correlated. 

In the noise-free case, £ K = £ T = 0, the path r is a perfect helix with radius r = «o/(^o + r o )> 
pitch 27ih = 2n t /(kq + Tq ), and helix angle 9 = tan~ 1 (/i /r ). We define the helix reference 
frame (R, hi, h 2 , h 3 ) by the linear transformation 



R = r + r n, h 3 = sin 6q t + cos 6*o b, 



(3) 



hi = — n and h 2 = h 3 x hi. Here, R(£) is the centerline of the helical path r(£) and h 3 is called the 
helix vector. The helix frame can be interpreted as the material frame of a solid disk with center R, 
see Fig. 1A: The disk translates and rotates such that a marker point on the disk's circumference 
traces the helical path r. For a perfect helix, R = Uh 3 , h 3 = 0, hi = u h 2 , h 2 = — c^ohi where 



2 



v = co h and u> = v (k% + t^) 1 / 2 is the frequency of helical swimming. The period of a helical 
turn is T = 2tt/uj . 

In the presence of fluctuations, R = vh 3 + t> ro(£ r b — £ K t). The helix vector h 3 performs a 
stochastic motion on the unit sphere which is characterized by 

(h 3 (0) • h 3 (t)} « exp(-t/t P ) (4) 

for times t longer than the correlation time of curvature and torsion fluctuations. In the following, 
we determine the persistence time tp in a limit of weak noise; the result is given in eqn. ©. The 
rotation matrix H(t) with Hki = hfc(O) • hi(t) is an element of SO(3). The Lie algebra of SO(3) 
is spanned by the infinitesimal rotations E - with (E -)jy = e k ji, j = 1, 2, 3. The time evolution 
of H(t) is given by a matrix-valued differential equation H = H • h with infinitesimal rotation 
h = cu E 3 + £jE - where we use Einstein summation convention for j = 1, 2, 3. From eqns. dHO), 
we find £ x = 0, £ 2 = ( r o£r — h ^ K ), and £ 3 = u (r £ K + h £ T ). The rotation of the helix frame 
after a time t consists of a rotation around h 3 (0) by an angle u t and random rotations around 
all axes due to the curvature and torsion fluctuations. We characterize these random rotations by 
continuous stochastic processes Ej(t) with 2, (0) = and write 

H(t) = exp(oo t Eg) • exp(H,E J ). (5) 

Note that exp(u; nT E 3 ) = 1 after n helical turns. The Ej represent generalized rotation angles: 
Si and S 2 describe rotations of h 3 . Symmetry implies (Si) = (S 2 ) = 0. We consider the limit of 
weak noise characterized by |5V(co>)| k, /v . One can develop a systematic expansion 

in powers of the noise strengths rj K , r\ T with r\ 2 K = v r dt\(€ K (0)£ K (t))\ and analogously rj T . 
We write = to denote equality to leading order in i] K , r\ T . For times t = nT longer than the 
correlation time of £ K , £ T but still shorter than T/(r) K + r) T ), we find 

(~j) = (~i) = 2Dt, (SiS 2 )^0 (6) 

with AD = S 2 (ujo) where 5 2 (o;) = uj 2 [hlS K (uj) + 7qSV(^) — 2r h Re S K>T ((v)] is the power 
spectrum of £ 2 . Hence, the stochastic motion of the helix vector h 3 can be effectively described as 
isotropic rotational diffusion with rotational diffusion coefficient D for long times. The derivation 
of © proceeds as follows: H(t) can be written as a time-ordered exponential integral H(t) = 
Texp f* dt'h(t'). To linear order in the noise strengths, S 2 + i S x = jjj dt' £ 2 (t) e^'^. Next, 
(El + S|) ^ J * dt x Si dt 2 (Uh) Uh)) e-^fr-fc) w 5 2 (^ ) Similarly, (S|) ^ S 3 (0) t where 
^(a;) is the power spectrum of £ 3 . Hence the swimming path r(t) is a noisy helix with a centerline 
R(£) that follows a persistent random walk (on time-scales larger than the correlation time of 



curvature and torsion fluctuations) Q2VU - This persistent random walk has a persistence time 



t P = (2D)- 1 = 2/S 2 (cu ) (7) 

that is governed by the power spectra of the curvature and torsion fluctuations evalutuated at the 
helix frequency cuo and a persistence length lp = vtp I119l l20n. 

A chemotactic chiral swimmer. We now consider a chiral swimmer in a concentration field 
c(x) of chemoattractant equipped with a feedback mechanism which allows it to dynamically 
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adjust its curvature and torsion in response to a chemotactic stimulus s(t). The stimulus s(t) = 
Y^j 8 {t—tj ) counts single chemoattractant molecules detected by the swimmer at times t j . The rate 
q = (s) of molecule detection of the swimmer is assumed proportional to the local chemoattractant 



concentration 112 lh 



(s(t))=q(t) = \c(r(t)). (8) 

When q is large compared to a typical measurement time a -1 of the swimmer and q{t) changes 
on a time-scale slow compared to the mean inter-event-interval 1/q, then we can replace s(t) by a 
coarse-grained version known as the diffusion limit 



s(t) « q(t) + VW)Ut) (9) 

where £ s (£) is Gaussian white noise with (£ s (£i)£ s (£2)) = <5(^i — In this limit f) < 1 where 
7) = (q(j)~ l l 2 characterizes the relative noise strength of s(t) for an averaging time a 112 ill . The 
chemotactic stimulus s(t) is transduced by a signaling system of the swimmer and triggers a 
chemotactic response which we characterize by a dimensionless output variable a(t) with a = 1 
for a time-independent stimulus s(t) = sq. We assume that aft) affects curvature and torsion in 



a linear way K,(t) = k q + Ki(a(£) — 1) and analogously for r tl8Q. Recall that swimming speed 



v{t) = v is assumed constant. For the signaling system relating stimulus s(t) and output aft), we 



use a simple dynamical system which exhibits adaptation and a relaxation dynamics fll8l l21. 



22] 



aa = ps — a, fip — p(l — a). (10) 

Here p(t) is a variable representing a dynamic sensitivity; a is a relaxation time and /i is a time- 
scale of adaptation. For a time-independent stimulus s(t) = sq, the system (flOl) reaches a station- 
ary state with a = 1, p — l/s . Small periodic variations of the stimulus s(t) = sq + s\ cosut 
evoke a periodic response of the output variable a(t) = 1 + s{Re Xa(^)e tuJt + 0(s\) with linear 
response coefficient Xa(u) = iu>fx/[s (l + iujfi — afiui 2 )]. 

Swimming in a concentration gradient. We consider a chemotactic chiral swimmer in a linear 
concentration field of chemoattractant 

c(x) = C + Ci ■ X. (11) 

Fig. IB shows an example of a stochastic swimming path r(t) in such a linear concentration 
field which has been obtained numerically. In the simulation, the chemotactic chiral swimmer 
detects individual chemoattractant moleculs arriving at random times (distributed according to an 



inhomogenous Poisson process with rate q(t) 112 ill ). 



We characterize the chemotaxis mechanism of a chiral swimmer in the limit where both 
chemoattractant concentration c is high with 77 = (Ac cr) -1 / 2 1, and the concentration gra- 
dient is weak with v = |ci|ro/co 1. The concentration gradient Ci is a sum ci = Cj|h 3 + 
with a component parallel to h 3 of length en = Ci ■ h 3 , and a component Cj_ = c x — C||h 3 perpen- 
dicular to h 3 of length c± = \c±\. While the swimmer moves in the concentration field along the 
noisy helical path r(£), the binding rate q(t) varies with time. In the limit of weak noise and a weak 
gradient, we approximate q(t) by the value for q(t) obtained for swimming along the unperturbed 
path with chemotactic feedback switched off q(t)X ~ c + cm (0)vt + c±(0)r cos(u t + (p ) where 
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ipo is the angle enclosed by cj_(0) and hi(0). It is this periodic modulation of q(t) which underlies 
navigation in a concentration gradient as it causes a bias in the orientational fluctuations of h 3 : 
The stimulus s(t) elicits a periodic modulation of the average curvature and torsion with ampli- 
tude proportional to c±. As a consequence, the expectation values (Si), (S 2 ) of the generalized 
rotation angles introduced in © are non-zero and scale with c±, (H 2 + iSi) = c± Ex a (uJo)e~ l{pQ t 
with e = XujqTo^qKx — r ri)/2. Similarly, (H 3 ) = c\\ It with e = fiu^ho^roKi + h Ti)/c UM- 



We can now derive an effective stochastic equation of motion for the helix frame in the limit 



rj, v <C 1 by a coarse-graining procedure as outlined in \lv\ . The Stratonovich stochastic differen- 
tial equation for the helix frame 

R =vh 3 

h 3 =-eRe[x a (w )c] + £ 1 h 2 -£ 2 hi (12) 
hj = -(hg • hj) h 3 + e kj3 uh k , j = 1, 2 



generates the statistics of the noisy helical path to leading order in v and rj 02811 - Here uJ = u + ec\\ 
and c = Cj_ + i h 3 x cj_. Eqn. (fT2l) contains a multiplicative noise term ^ x h 2 — £ 2 h! where ^ 
denotes Gaussian white noise with (^k{ti) Ciih)} = 2D5ki5(ti — t 2 ). Here D plays the role 
of a rotational diffusion coefficient and is given by D = \e Xa(^o) / r o\ 2 Co/A. Note that -D is 
concentration dependent with D ~ 1/cq. In the deterministic limit = = 0> we recover the 



results from [18]. Eqn. (fT2l) provides a coarse-grained description of the time evolution of the 
helix frame on time-scales larger than the correlation time a of curvature and torsion fluctuations. 

Effective dynamics of the alignment angle. In a linear concentration field ([111) . the quantity of 
interest is the alignment angle ip enclosed by the helix vector h 3 and the direction of the gradient 



Ci 118D, see Fig. 1A. The symmetries of the problem imply that the dynamics of ip decouples from 
the other degrees of freedom of the helix frame. From (fT2)) . we find by using the rules of stochastic 
calculus 

ip = -/3 sin + £ + cot (13) 

Here £ denotes Gaussian white noise with (£(£i)£(t 2 )) = 2D 5{t\ — t 2 ). The alignment rate 
is P — — |Vc| e ReXa(^o)- In me absence of fluctuations, we recover the deterministic limit 



ip = —(3 smip Ill8h . In this limit, the steady state is characterized by either parallel alignment of 
helix vector and concentration gradient with ip = for j3 > or by anti-parallel alignment with 
ip = Ti for (3 < 0. Eqn. (TT3l) contains a noise-induced drift term D cot ip which diverges for ip = 
and ijj = ir implying that noise impedes perfect parallel or anti-parallel alignment of the helix 
vector. 

The corresponding Fokker-Planck equation for the probability distribution P(z, t) of z = cos tp 
with \z\ < 1 reads P = —d z [(l — z 2 )((3 — Dd z )]P. Fig. 1C compares P(z,t) to a histogram 
of z obtained from simulating 10 5 chemotactic chiral swimmers in a linear concentration field. 
The distribution P(z,t) relaxes to a steady state distribution Pq(z) ~ exp(/3z/ D) on a time- 
scale which is set by the inverse alignment rate fi^ 1 . This steady-state distribution Pq(z) has its 
maximum at z* = ±1 for j3 ^ 0, respectively. The first moment of Pq{z) is given by the Langevin 



function |23|] 



(z) = coth(Pe) - Pe" 1 , Pe = f3/D (14) 

where Pe describes a Peclet number of rotational motion. Note that this result for the mean orienta- 
tion of a chemotactic chiral swimmer is formally equivalent to the orientation of a polar molecule 
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FIG. 1 : A. A helical path r can be described as the trajectory of a point on the circumference of a disk which 
rotates and translates; see text for details. B. Helical swimming path r of a chemotactic chiral swimmer 
in a linear concentration field: The helix vector b.3 fluctuates around the direction of the concentration 
gradient Vc. Parameters were vq = tq/g, tq = 0.2/tq, —k\ = t\ = 0.5/ro, /U = a, Xcq = 10 2 /a, 
|ci| = 10~ 2 co/?"o- C. Histograms of z = cosip where tp equals the angle enclosed by the helix axis and 
the gradient direction for a simulated ensemble of helical swimming paths as in B with initial distribution 
P(z, 0) = 5{z) at times t = 10T, 100T. Also shown is the analytical solution P(z, t) (red). Approximately, 
P(z, 100T) equals the steady-state distribution Pq{z) ~ exp(Pez). 



in an external electrical field: Eqn. (fl~4l) with Pe replaced by |m| \Fi\/(k B T) also describes the 
mean orientation (z) = (m ■ E)/(|m| |E|) of a polar molecule with dipole moment m subject to 
rotational Brownian motion in an electric field E ll23ll . Note that eqn. (fl~4l) characterizes an active 
process while a polar molecule is an equilibrium system. 

At steady state, a chemotactic chiral swimmer moves up a concentration gradient with average 
speed (z)v. The chemotaxis index CI is defined as the ratio of this average speed gradient-upwards 
and the swimming speed v 



CI = (z)CI r 



CL 



v/vq = sin Q . 



(15) 



Note that CI approaches its maximal value CI max for Pe > 1. This condition is satisfied already 
beyond moderate concentration gradients with | Vc| > |£|/(Ar ) 2 . The maximal value CI max for 
the chemotaxis index is limited only by the geometry of helical swimming. 

Relation to experiments. Chemotaxis of sperm cells has been extensively studied for sea urchin 
sperm cells tl6fl. Tracking experiments in three dimensions show that these sperm cells swim 
along noisy helical paths with typical values for swimming speed, average curvature and torsion 
v w 100-200/ims" 1 , k w 0.025-0.05 /im -1 , r w -0.0025 /im" 1 For comparision, the 



length of the sperm tail is L ps 50 jttm [|15|1 . Using a two-dimensional experimental setup in which 
sperm cells swim along a circular path, it has been shown that a periodic chemotactic stimulus 
causes a phase-locked periodic swimming response [|15L \24\ . Such a behavioural response is 
consistent with our model of a chemotactic chiral swimmer. 

In a pioneering experiment, C. J. Brokaw observed helical swimming paths of bracken fern 
sperm cells in a shallow observation chamber [0,3] [ 29 ] . In the absence of chemoattractant, sperm 
swimming paths were noisy helices whose centerlines could be described as planar persistent 
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random walks with persistence time tp^d ~ 5s and net speed v ~ 200 /ims^ 1 , corresponding 
to a persistence length of Ip^d ~ 1 mm - Accordingly, the planar orientational fluctuations of the 
helix vector are characterized by a rotational diffusion coefficient D = t~p l 2d = 0.2 s" 1 . In a 
strong concentration gradient of chemoattractant, sperm swimming paths were bent helices which 
aligned with the gradient direction at a rate proportional to the relative strength of the concentration 
gradient (3 ~ 150 /zms -1 | Vc|/c yOfl . In an initially homogeneous concentration field of charged 
chemoattractant, alignment of helical sperm swimming paths could also induced by applying an 
external electrical field |E|. In this case, it was found that the alignment rate is proportional 
to the field strength (3/\E\ ps 1.6s" 1 (V/cm) -1 . The mean alignment (z) = (h 3 • E)/|E| of 
helical paths at steady state was measured as a function of field strength |E|. The experimental 
data could be well fitted by eqn. ([T4|) assuming Pe ~ |E| and yielded Pe/|E| ~ 8 (V/cm)~ 
Olll . The above estimates for D and /3/|E| give approximately the same value for Pe/|E| = 
{(3/D)/\Ei\ The physical origin of helix alignment in an electrical field is not entirely known: 
The observed alignment might be due to electrohydrodynamic effects resulting from sperm cells 
binding chemoattractant ions (with sperm cells effectively behaving as electric dipoles) [6i]. An 
alternative possibility is that the electric field induces a concentration gradient of chemoattractant 
ions and that the observed alignment of helical paths is a result of chemotactic navigation in this 
gradient. 

Conclusion. In this Letter, we studied the stochastic differential geometry of noisy helical 
swimmingpaths which is relevant for many biological mircoswimmers with chiral propulsion 
[flBSHlEli, 10, U, 12]. A simple feedback mechanism enables a chiral swimmer to navigate 
along a helical path upwards a concentration gradient of chemoattractant. Chemotaxis alongnoisy 
helices is employed by sperm cells and possbily other biological microswimmers lalal?],lEl 9 l ll2|l . 
A similar mechanism underlies phototaxis of the unicellular flagellate Chlamydomonas P5Q. and 
is found in phototactic marine zooplankton [| 1 0L 1 1 lf l - Our theory shows that navigation along he- 
lical paths is remarkably robust in the presence of fluctuations: An effective rotation of the helix 
vector is determined by integrating its orientational fluctuations over several helical turns. Con- 
sequently, a small bias in these orientational fluctuations due to chemotactic signaling results in 
robust steering and the helix vector tends to align with the concentration gradient Vc. If chemo- 
tactic signaling is adaptive, the alignment rate (3 is proportional to the relative strength of the 
concentration gradient | Vc|/c. After a transient period of alignment of duration a chemotac- 
tic chiral swimmer moves upwards the concentration gradient with an average speed that is only 
limited by the geometry of helical swimming provided the strength of the concentration gradient 
exceeds a characteristic value. We conclude that temporal sampling of a concentration field along 
a helical path provides a universal strategy for chemotaxis which is highly adapted for a noisy 
environment. 
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Eqn. (fT2l generates stochastic processes as denned in ((5]) whose second-order statistics is correct to 
leading order in v and r\. 

Bracken fern sperm cells are very different from animal sperm cells, but also motile with thrust gen- 
erated by several fla gella . 



The value quoted in H 1 811 was incorrect due to a confusion of decadic and natural logarithm. 
Swimming of sperm cells was restricted to a shallow observation chamber of height 100 /jm which 
affects the statistics of z: If the helix vector is constrained to a plane parallel to E, our theory predicts 



(z)2d = ^i(P e )/-^o(Pe). This also describes the data in (60 and gives Pe/|E| « 5 (V/cm 
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